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The double beta decay of ^^^Mo to the ground state and excited states of 
^ ' ^^^Ru is analysed in the context of the pseudo SU(3) scheme. The results of 

this deformed limit are compared with the vibrational one based on the QRPA 
formalism. Consistency between the deformed limit and the experimental in- 
0^ ' formation is found for various /5/5 transitions, although, in this approximation 



some energies and B(E2) intensities cannot reproduced. 
PACS numbers: 21.60.Fw; 24.80.Ba; 27.60. +j 

The two neutrino mode of the double beta decay {f3f32u) is an exotic but 
allowed second order process in the standard model. It has been detected in 
nine nuclei and has served as a tool to develop the most refined experimental 
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detection methods for very low signals. It also represents a severe test on the 
nuclear structure description of those nuclei j^ . 

The ^^^Mo is an excellent subject of study. Its large Q-value , Q^^ = 
3.034:MeV, favors its decay. Four groups have reported on the detection of 
the P(32u of ^^^Mo to the ground state of ^^^Ru with half-lives of the order of 
lO^^yrs ||2[^,QJ5| . Also the decay to the first excited 0^ state has been reported 

©I- 

This abundance of experimental information allowed a careful check of 
the nuclear models used to describe the two neutrino mode. The QRPA 
model, both in its conventional and number projected versions, was found 
to collapse when realistic particle-particle interactions are used [0,^]. Under 
the assumption that the excited 1130 keV 0^ state in ^^^Ru is a member 
of a 0^ , 2+ , 4+ two-phonon vibrational triplet the QRPA calculation exhibits 
an overestimation of the amplitude of the P(32u decay to this excited state. 
However the decay to the ground state is reasonably well reproduced by using 
the same approximation. The single beta-decay transition between the initial 
and intermediate nuclei is also overestimated by the theory |p. Assuming 
realistic nucleon-nucleon interactions similar results were found, i.e. the Pf32u 
decay amplitudes are overestimated. This result is particularly true for the 



decay to the first 2~^ state [10| 



All the works dealing with the /?/32!/ decay of ^^^Mo found a strong dom- 
inance of one single particle transition: {97/2(1^'))'^ — > (fi'9/2(p))^- This lack of 
collectivity was mentioned as a possible cause of the failure of the QRPA to 
describe the P(32u decay of this specific nucleus |^] while having success in 
others 0Jl3,|l2|,|l3|. 



Other complementary explanation could be that deformation effects were 
destroying the one-phonon-two-phonon structure in ^^^Ru, and substituting it 
with rotational bands built upon a deformed ground state or some vibrational 



intrinsic excitations [|T^]. In the present paper we explore both possibihties 
using the pseudo SU(3) formahsm. 

The transition from spherical to deformed shapes in Mo isotopes was 
study microscopically several years ago, commenting on the deformed char- 



acter of some first excited 0"*" states [14]. In a similar way within the General 



Collective Model (GCM) [^ a study of the potential energy surfaces (PES) of 
the Ru isotopes |1^ exhibited a transition from a harmonic oscillator (^^Ru) 
to a nucleus with triaxial dominance i^^^Ru). All Ru isotopes between ^^^Ru 



and ^^'^Ru behave very similar |16|. The PES has a global minimum at /? = 0, 
suggesting an harmonic oscillator for the low lying states. There is also a 
local minimum at /3 approximately 0.4 and 7 = 24", which suggests triaxial 
shapes for excited states. However the energy of the ground state lies above 
the potential barrier between the spherical and the triaxial minimum. As 
a consequence, the ground state, the first and second excited 2"*" state are 
spread out between the global and triaxial minimum showing the behavior 
of an anharmonic oscillator, while the first excited 0^ and third excited 2+ 
state are peaked within the triaxial minimum. As a picture we obtain the 
spectrum of an anharmonic oscillator were the first excited 0"^ state does not 
fit into it. Nevertheless the B{E2) transitions follow the characteristic of an 
anharmonic oscillator within a simple consideration, i.e., the transition from 
the first exited 0^ state to the first excited 2+ state is strong compared to 
the second excited 2"*" state, which is weak. This did lead in the past to the 
assumption that "^^^Ru resembles more an anharmonic oscillator. 

We will use the pseudo SU(3) model |17] to provide a microscopical de- 
scription of the ground and excited low energy states in ^^^Ru, and the ground 
state of ^^^Mo including deformation. This exercise will supply a rotational 
limit in the description of these nuclei, to be compared with the spherical, 
vibrational limit obtained with the QRPA. 



The double beta decay, when described in the pseudo SU(3) scheme, is 
strongly dependent on the occupation numbers for protons and neutrons in 
the normal and abnormal parity states n^, n^, n^, n^ [18|. In particular, the 



Pf32u decays are allowed only if they fulfil the following relations: 



</ = <. + 2, </ = <. </ = <.' </ = <-2 (1) 

These numbers are determined by filling the Nilsson levels from below, as 



discussed in [18|. They are exhibited in Table 1 as the deformed occupation 
(def-o) case. Both ^^^Mo and ^^^Ru have deformations P ~ 0.23. This 
number must be taken carefully, given these nuclei do not exhibit a rotational 
spectra, but is a hint about their average deformations. We were forced to 
select a slightly higher deformation for neutrons than for protons. If not, 
the removed pair of neutrons in ^^^Ru would come from the abnormal parity 
orbital, and in this case the expressions (|I|) would not be fulfilled suppressing 
the Pf32u decay. 

In the third and fourth row of Table 1 the occupations in the spherical 
occupation (sph-o) limit (/? < 0.09) are given. This is a very interesting 
limiting case, in which all the proton holes belong to the abnormal gg/2 orbit 
and all the neutron particles to the normal parity orbitals. 

In the abnormal parity space only seniority zero configurations are taken 
into account. This is a very strong assumption, which in future works is ex- 
pected to be improved, but is quite useful in order to simplify the calculations. 
Its effects in the present calculation are discussed below. 

Under the assumption of a Hamiltonian consisting in a Nilsson mean field 
and a strong quadrupole-quadrupole interaction, the deformed wave functions 
of ^^^AIo and ^^^Ru, with angular momentum J and belonging to the ground 
state {a = g.s.) or the excited (a = exc.) band can be constructed according 
to 111. The result is: 



\^o^Mo{Ru),J+a)= |{2^}(A^,^^); {2^}(A,,/x,); l(A^),i^ = IJ >iv 
I (59/2)"% J^ = 0; {hn/2r^, J^ = 0>A , 

Each (X'sjfi's) are also depicted in Table 1. In this approach, instead 
of assuming the first excited 0'*' and the second 2"^ states as parts of a two 
phonon triplet, together with the first 4"'" states, it is assumed that the first 
0^, 2^, 4^ states of ^^^Ru are the low energy sector of a rotational band based 
in the normal (A, fJ-)g.s. = (8, 6) strong coupled pseudo SU(3) configuration 
and that the second 0"*", 2+ excited states belong to a second rotational band 
described by the (A,/i)exc. = (10,2) pseudo SU(3) configuration. Note that 
in the sph-o case the normal parity orbitals (A^ = 3) are filled, giving rise 
to the irrep (0,0). It implies that these protons are dynamically inert. It 
also has the consequence that there exists only one strong coupled stated 
(0, 0)7r X (12, 0)1, = (12, 0), and in this approximation there is not a 0"*" excited 
state in our restrictive Hilbert space. 

In order to analyse the spectra and transitions amplitudes of ^^^Ru we 



(2) 



have selected a simplified version of the pseudo SU(3) Hamiltonian |17], i.e. 



/7 = ^F,-ixQ"-Q" + aL2 a = 7r,u, (3) 

a 

where the algebraic quadrupole operator Q** = J2s il-rrs + Qfs} ^-^ts only 
within a shell and do not mix different shells and 

where rja = fja -\- 1 denotes the phonon number operator, fito determines the 
size of the shell and a = t: ox u. A constant term \y (V,r) is included, which 
represents the depth of the neutron (proton) potential well. 

The abnormal parity sector of the wave function contributes with a con- 
stant to the eigenvalues of this Hamiltonian, because the seniority zero approx- 
imation inhibits any dependence with the total angular momentum. Thus, the 



expectation value of the Hamiltonian between the ^^^Ru wave functions only 
depends on the (A, /i) and J = L of the normal parity orbitals. Explicitly: 

<i°o Ru, (A, fi)K = 1, L = J M\H\^^^Ru, (A, fi)K = 1,L = JM>= 
N.hio - 2xC2 - (|x + a)L{L + 1) 

In the above expression Ng represents the number of proton plus neutron 
phonons, which depends only on the normal parity occupations n^ and n^ 
and do not depend on the irrep or on the angular momentum L. The eigenval- 
ues of the second order Casimir operator of SU(3) isC2 = (A+/i+3)(A+/u)— A/u. 
The irrep (8, 6) has the maximum value of C2 = 190 and yields the ground 
state band. The second 0+ and 2"*" states were associated with the irrep 
(10,2), withC2 = 160. 

The two parameters of the Hamiltonian (3) are fitted by means of the 
excitation energies of the first 2"'" and the first excited 0"*" states. The results 
are x = 18.83keV, a = 61.75keV, which are reasonable numbers as an inter- 
polation between those needed to reproduce the rotational spectra of '^'^Mg 
|p!9| and ^38^/ ^2C\. In the comparison it must be taken into account that the 



Refs. |19| and [|2^ use symplectic extensions of SU(3) model and thus they 
need smaller coupling constants. 

In Table 2 the experimental and theoretical low lying states in ^^^Ru are 
presented, with the angular momentum, parity and energy stated in each 
level, as well as the corresponding irrep (A, fi) of each band. As expected, the 
energies of the first 4^ and the second 2^ states are overestimated. This is a 
remainder that the present exercise is a rotational limit of a triaxial nucleus, 
whose low energy spectra evokes an anharmonic oscillator, as it was mentioned 
above. 

The average quadrupolar moment Qo is related with a particular B(E2) 
value as follows: 



(5) 



S(£;2,0+^2+) = ^|Qo|' (6) 

In the pseudo SU(3) scheme, this intrinsic quadrupolar moment is easily 
evaluated as ||2^ 

Qo = el"Q^ + elffQ, Q„ = ^^i±l(2A„ + /x„) (7) 

where Qq is given in units of e 6q with 6q = -p^ = 4.70/m^. Using the wave 
functions of the ground state band (^) the intrinsic quadrupole moments are: 

Q^ = 5.33 Q^ = 22.5 (8) 

The experimental value is Qo{exp) = 57. 7e b^. To reproduce this result 
unreasonable effective charges e^^^ « 2.9e and e^-^-^ ^ 1.9e would be needed. 
This failure of the pseudo SU(3) model in reproducing the quadrupole mo- 
ment of ^^^Ru is related with the use of the seniority zero approximation 
for nucleons in abnormal parity orbitals. In the case of ^^^Ru case even 
in the def-o limit there are only two normal holes and the associated irrep 
(0, 4) gives a very little proton quadrupole moment. The greater part of the 
proton quadrupole moment must come from protons in the abnormal parity 
orbital gg/2, whose contribution is neglected in the seniority zero approxima- 
tion. These limitations were discussed previously [^,22]. In this picture large 



effective charges are always needed in order to compensate this effect. For 
the case of ^^^Ru these effective charges are too large. Although we have pre- 
sented here some limitations of the rotational model of ^^^Ru, it is successful 
in describing the /?/?2i/ decay of ^^'^Mo. 

We turn now to the study of the two neutrino mode of the double beta 
decay /3/32i/ of ^^^AIo into the ground state, the first 2^ and the first excited 
0"^ states of ^^^Ru. The mathematical expressions needed to evaluate the 



/3/?2!/ to the ground state of Ru can be found in |18]. The same formulae 



works for the decay to the first excited 0"^ state, replacing the strong coupled 



irrep (8,6) by (10,2). The decay to the first 2+ state requires a different 
expression, which can be found in [^. The formulae for this decay resembles 
that of the decay to the 0^ states, but the energy denominator is elevated to 
the third power. Being in general this energy of the order of 10 AleV this 
power implies a factor 100 of suppression for this matrix element JIJ^lOl. 



In Table 3 the dimensionless /?/?2i/ matrix elements for the decay of 
lOOjy^Q to the ground state, the first 2+ and the first excited 0"*" states of 
^^^Ru are presented. The first two rows include the results of the QRPA cal- 
culations obtained using a (5-interaction y] and an interaction derived from 
the Bonn one-boson-exchange potential using G-matrix techniques |l^. The 
third row shows the results for the pseudo SU(3) calculation using the def-o 
wave functions given in Eq. (|2|) . The energy denominator used for the decay to 
the ground state, determined by fixing the energy of the Isobaric Analog State 
\l^is E = Eq- 2hujk^{r]^ + I) + Ac = 7.95MeV, where Eq = ^ Qpp + rUeC^ 



and Ac is the Coulomb displacement energy. The matrix elements are given 
in units of the first or third power of the inverse electron mass Q , since the 
energy denominators have been divided by the electron rest mass. For the 
first 2"*" state the energy denominator is equal to l.QSMeV and for the first 
excited 0"^ it is 7.39MeF. The results for the sph-o pseudo SU(3) approach 
are given in the fourth row. As it was anticipated in the spherical occupation 
limit the decay to the excited 0"*" state is absent. 

The experimental matrix elements showed in Table 3 were extracted from 

1/2 

the measured half-lives T2I using phase space integrals Ggt^ which were 
obtained following the prescriptions given in |24| with qaIsv = 1-0. These 
data are reproduced in the last two rows of Table 3. 

It is remarkable that the sph-o pseudo SU(3) limit reproduces quite well 
the experimental matrix element for the (3(32y to the ground state and that 
the def-o pseudo SU(3) limit is able to do the same with the matrix element 

8 



related with the (3(32u decay to the first 0+ state. Both the abihty to reproduce 
the numbers and the difficulty to fit both in the same model are pointing to 
the description or ^^^Ru as a triaxial nucleus. 

The P(32u matrix element for the 0^ — > 2+ decay is strongly cancelled in 
comparison with the QRPA one [0]. In the sph-o case, which successfully 
reproduce the 0+ — > 0^{g.s.) decay, the reduction has two contributions: i) the 
matrix element of the operator ui • (T2 which connects the 0"'' states is a factor 
ten greater than the matrix element of the operator [ui 'S>cr2]'^ which mediates 
the other decay; ii) the energy denominator in the pseudo SU(3) model is 
approximately three times greater than the energy denominator associated 
with the first intermediate state 1^, giving a 3^ = 27 additional reduction 
when compared with the QRPA result. 

In order to estimate the effect of the quadrupole-quadrupole correlations 
over the Pf32u matrix elements, we have performed also a simple estimation 
of these matrix elements under the pure seniority zero approximation. It is a 
model in which only the gg,2 and §^,2 orbitals are active, with the nucleons 
paired to angular momentum zero. There are 2 protons and 8 neutrons in 
^^^Mo and 4 protons and 6 neutrons in ^^^Ru. These two (four) protons are 
described in the same way in this simple approach and in the sph-o pseudo 
SU(3) scheme but neutrons have strong mixing with their partners in the 
A^ = 4 oscillator shell in the pseudo SU(3) scheme which is absent in the 
other. 

The seniority zero matrix elements [|l^] are exactly the same as the pure 
pairing ones (7,13|, i.e. 

M2M - 0+(5...)) = E ^^^"^"^P < ^";"^ >' (9) 

p„ ep + en- i/o 

There exist only one (p, n) = (59/21 57/2) configuration in this approach, 
limiting the sum to only one term. The pairing occupation numbers w^ in 



the initial (unbarred) and final (barred) nuclei are u^ = 0.8, v^ = 0.4, f^ = 
1.0, n^ = 0.25. Using this numbers together with the energy denominator 
given above we have obtained the matrix element for the P(32u in this seniority 
zero approach which is given in the fifth row of Table 3. The comparison with 
the pseudo SU(3) matrix elements exhibits a reduction by a factor two in the 
sph-o limit and by a factor three in the def-o one. This reduction is due to the 
quadrupole-quadrupole correlations, which strongly mix the §^,2 with their 
shell partners, and also mix these neutrons with the normal parity protons. 
Using an energy denominator consisting essentially of twice the pairing gap 
plus half the Q-value (~ 3MeV) would increase the seniority zero matrix 
element making it worst when compared with the experimental value. 

In the present paper we have studied the Pf32i/ of ^^^Mo to ^^^Ru in the 
context of the pseudo SU(3) model. It was also mentioned that in the geo- 
metrical model, while the ^^^Ru ground state essentially feels an anharmonic 
oscillator potential, the first excited O'^ state is concentrated in the deformed 
region. We have explored the possibility of a coexistence of shapes, trying to 
find out a possible explanation for the difficulties found in previous attempts 
to describe the (3(32u of ^^^Mo to excited states of ^^^Ru which described the 
first excited O"*" state as a two phonon state. 

We have used the pseudo SU(3) model as a deformed limit and compared 
it with the QRPA, which is taken as the spherical limit. Under this scheme, 
and with the additional approximation of a seniority zero wave function for 
the abnormal parity nucleons, we have generated a rotational spectrum for 
^^^Ru which does not resemble so much the experimental one. Also, we 
were not able to generate enough proton quadrupole moment since we have 
used the seniority zero approximation for abnormal parity nucleons. Besides 
this limitations, the description of the Pf32u was successful. We were able to 
reproduce the experimental nuclear matrix elements, without any parameter, 
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but using the deformed occupation pseudo SU(3) wave function for P(32i, to 
the first excited 0"*" state, and the spherical occupation pseudo SU(3) wave 
function for the f3f32u to the ground state. Fitting both decays in the same 
model is not possible, at least within the present scenario. 

We acknowledge P. Vogel for pointing out an incorrect rescaling in the 
phase space integrals and S. Pittel for interesting comments about the micro- 
scopic origin of deformation in ^^^Mo and about the seniority zero approach. 
This work was supported in part by Conacyt under contract 3513-E9310, and 
by a Conacyt-CONICET agreement under the project Double Beta Decay. 
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Table captions 



Table 1 Occupation numbers and (X's, fi's) for the deformed occupa- 
tion (def-o) and spherical occupation (sph-o) pseudo SU(3) wave function for 
^^^Mo and ^^^Ru. 

Table 2 Experimental and theoretically determined energies (in keV) for 
the ground state (g.s) end excited (exc.) band of ^^^Ru. 



Table 3 The dimensionless Pf32u matrix elements for the decay of 
^^^Mo to the ground state, the first 2^ and the first excited 0"^ states of 
^^^Ru, evaluated with different models. The experimental matrix elements, 
measured half lifes T2J and phase space integrals G'-^^ are also listed. 
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Table 1 







ni n^ 


N N 


(A7r,M7r) 


{X,y,l^u) (A,^)g.s. 


(A, ^)ezc. 




Mo 


4 2 


10 6 


(0,4) 


(12,0) (12,4) 




def- 


- 














Ru 


6 2 


10 4 


(0,4) 


(8,2) (8,6) 


(10,2) 




Mo 


2 


12 8 


(0,0) 


(10,4) (10,4) 




sph - 


- o 














Ru 


4 


12 6 


(0,0) 


(12,0) (12,0) 





Table 2 





5.5. 


band 


exc. hand 


J+ 


th. 


exp. 


th. exp 


irrep 


(8,6) 




(10,2) 


0+ 








1130 1130 


2+ 


540 


540 


1362 1670 


4+ 


1800 


1227 
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Table 3 





0+ ^0+(5.s.) 


+ ^0+(exc.) 


0+^2+ 


QRPAg 




-0.256 


-0.256 




QRPA |1C1 




0.197 


-0.271 


-0.033 


pseudoSU {3){def - 


-o) 


-0.108 


0.098 


1.53 X 10"^ 


pseudoSU {2,){sph - 


-o) 


0.152 




7.3 X 10-5 


seniority 




-0.323 






experiment 


±0.150 


±0.092 


< 0.106 


Tli' [lO^^r] 




1.15 ii 


178 |] 


> 11.5 [^ 


Ggt [10-2%r-i] 




387 


6.61 


76.6 
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